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ABSTRACT 

Context. The Deep Extragalactic VLBI-Optical Survey (DEVOS) aims at constructing a large sample of compact radio 
sources up to two orders of magnitude fainter than those studied in other Very Long Baseline Interferometry (VLBI) 
surveys. Optical identification of the objects is ensured by selecting them from the Sloan Digital Sky Survey (SDSS) 
list. 

Aims. While continuing to build up the DEVOS data base, we investigated how the VLBI detection rate could be 
enhanced by refining the initial selection criteria introduced in the first paper of this series. 

Methods. We observed 26 sources in two adjacent, slightly overlapping 2°-radius fields with the European VLBI Network 
(EVN) at 5 GHz frequency on 2 March 2007. The phase-reference calibrator quasars were J1616+3621 and J1623+3909. 
The objects selected were unresolved both in the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey 
catalogue and the SDSS Data Release 4. 

Results. We present images of milli-arcsecond (mas) scale radio structures and accurate coordinates of 24 extragalactic 
sources. Most of them have never been imaged with VLBI. Twenty-two compact radio sources (85% of our initial 
sample) are considered as VLBI detections of the corresponding optical quasars in SDSS. We found an efficient way to 
identify quasars as potential VLBI targets with mas-scale compact radio stucture at > 1 mjy level, based only on the 
FIRST and SDSS catalogue data by applying simple selection criteria. 

Key words, techniques: interferometric - radio continuum: galaxies - galaxies: active - quasars: general - surveys 

1. Introduction for more than 3000 unresolved radio sources from the Faint 

Images of the Radio Sky at Twenty-centimeters (FIRST) 

53 . Distant extragalactic radio sources detected with Very Long gurvey catalogue ( White et aL \m7p\ with integrated flux 

Baseline Interferometry (VLBI) are typically active galac- density of at least 20 mJy xhege sources could be targeted 

tic nuclei (AGNs): the rest-frame brightness temperature with phase-referenced radio interferometric observations by 

of a 1-mJy source radiating from a region subtending a few selecting nearby bright and com p a ct calibrator sources from 

milli-arcseconds (mas) at 5 GHz is at least 10 b - 10' K, of- e g the yLBA Calibrator Survey (VCSjE To indicate the 

ten much higher. The aim of the Deep Extragalactic VLBI- pot ential of a FIRST- and SDSS-based VLBI imaging sur- 

Optical Survey (DEVOS; Mosoni et al. WM Paper I here- vey we note that optical counter p ar ts of at least 30% of all 

after) is to build up a large sample of optically identified FIRST S0U rces are expected to be found in SDSS (Ivezic 

AGNs imaged with VLBI at mas angular resolution. Many et al [^jrj^ A recent sensitive single-baseline Arecibo- 

of the potential astrophysical and cosmological utilisations Effelsberg V LBI survey of ~ 1000 FIRST sources indicated 

of such a sample require the optical identification and spec- that about one third of all sources had mJy-level compact 

troscopic (or photometric) redshift of the sources. These stucture; irrespective of the total flux density (Porcas et al. 

data can primarily be obtained from large optical sky sur- |2QQ4[) 

veys. In particular, the latest released] of the Sl oan D igital The technical feasibility of DEVOS has recently been 

Sky Survey (SDSS; Adelman-McCarthy et al. |2007J con- demonstrated in a pilot experiment (Paper I). However, at 

tarns imaging data for over 20 000, and spectroscopic data present reac hi n g the ambitious goal of increasing the total 
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number of VLBI-imaged extragalactic sources to 
Send offprint requests to: S. Frey very demanding in terms of observing resources. This led 
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us to somewhat modify the sample definition criteria used 
in Paper I in order to increase the throughput of the survey. 
We studied two adjacent celestial fields, each of 2° radius, 
to search for mas-scale structures of radio-emitting quasars 
detected with the SDSS. Here we describe the target se- 
lection (Sect. [2]), the phase-referenced VLBI observations 
(Sect. [3]) and the data reduction procedure (Sect. QJ. The 
results are presented in Sect. [5] The VLBI images of the 
target sources are available in the online version of this pa- 
per. Implications of our results are discussed in Sect. [51 

2. Sample selection 

The survey approach for the DEVOS pilot project (Paper 
I) was to select unresolved sources from the FIRST sur- 
vey (White et al. 11997ft with integrated flux density Si. 4 > 
30 mJy in the vicinity (within 2°) of a selected cali- 
brator source. Then 5-GHz phase-referencing observations 
with the UK Multi-Element Radio Linked Interferometer 
Network (MERLIN) were performed to pinpoint potential 
compact sources suitable as VLBI imaging targets. In the 
end, 19 sources (40% of the original FIRST sample) could 
be detected with a global network of VLBI antennas at 
5 GHz. These DEVOS pilot observations showed, among 
others things, that all the 6 FIRST sources that appeared 
stellar in SDSS imaging (i.e. optical quasars) could be de- 
tected and imaged with VLBI. In the case of these objects, 
the offsets between the VLBI and SDSS positions were 
between 30 and 90 mas. Independent comparisons of ac- 
curate radio interferometric source positions with the op- 
tical c oordina tes of their SDSS counterparts (Frey et al. 
I2006al[2006bl Lambert et al. l2lTOrJj) also indicate that SDSS 
coordinates of radio-emitting quasars are accurate enough 
to be used as a priori values for phase-referenced VLBI ob- 
servations. 

Building on this experience, we modified the sample se- 
lection criteria applied in Paper I, concentrating on SDSS 
quasars only. The reason is to find a way to speed up the 
construction of the radio-optical data base by means of (i) 
reducing the radio observing resources required, in partic- 
ular eliminating the need of "filtering" MERLIN observa- 
tions, and (ii) considerably improving the detection rate of 
the VLBI observations. These improvements are achieved 
by tightening the original selection criteria. This comes to 
an expense of losing some of the generality of the DEVOS 
approach and avoiding potentially interesting individual 
objects or classes. We stress, however, that the survey can 
be supplemented with the sources missed by these more 
stringent criteria at a later stage. 

With this compromise, the objectives of many poten- 
tial scientific applications of DEVOS outlined in Paper I 
can still be addressed. In particular, the morphology of 
low-luminosity quasars could be compared with that of the 
high-luminosity ones; cosmological tests could be conducted 
with substantially larger samples; precise astrometric posi- 
tions could be used to directly link the radio and future 
optical reference frames; the sources detected could serve 
as phase-reference calibrators - either in the traditional 
"nodding" style or in the in-beam mode - for other VLBI 
projects. 

The target extragalactic radio sources for the observa- 
tions described here were chosen from the FIRST survey 
data base obtained with the US National Radio Astronomy 
Observatory (NRAO) Very Large Array (VLA) at 1.4 GHz 



(White et al. I1997|) . Their optical identification is ensured 
by cross-comparing the list with the 4th Data Release 
(DR4 of the SDSS (Adelman-McCarthy et al. [2006)1 - the 
actual release at the time of preparing for the observations. 
The selection criteria were as follows: 

1. the integrated flux density at 1.4 GHz is Si. 4 > 20 mJy; 

2. the angular size of the source is 8 < 5" (i.e. unresolved 
with the VLA in the FIRST); 

3. the angular separation of the target source from the 
phase-reference calibrator source selected is less than 
2°; 

4. the source has an unresolved (stellar) optical counter- 
part in SDSS DR4 within l'/5 from the FIRST position. 

Note that the minimum flux density is lower here than 
that applied in Paper I (30 mJy). Very few Galactic radio 
stars would satisfy our selection criteria, thus the extra- 
galactic nature of the selected objects is practically ensured 
and could be verified with SDSS spectral data. 

Radio interferometric images of compact sources that 
are up to 2 orders of magnitude weaker than currently avail- 
able in other large-scale VLBI surveys can be made with 
the technique of phase-referencing (e.g. Lestrade et al. 1990; 
Beasley & Conway Il995p . It is based on extending the co- 
herence time by using regularly interleaving observations 
of a nearby bright and compact calibrator source. Here we 
select suitable calibrators first, then look for the potential 
target objects that are found in their close vicinity in the 
sky. Thus single calibrators can serve as phase-reference 
sources for several targets. This method was introduced by 
Garrett & Garrington ([1998]) and Garrington et al. (fl"999"l) . 

We choose two bright compact radio sources 
(J1616+3621 and J1623+3909) as phase-reference cali- 
brators. These are known VLBI sources (Beasley et al. 
120021 Helmboldt et al. 12007ft and have a relatively large 
number (13 and 16) of matching FIRST/SDSS quasars in 
their neighbourhood, according to the selection criteria 
above. The equatorial coordinates of the calibrators are 
determined in the International Celestial Reference Frame 
(ICRF) with mas or sub- mas accuracy (Petrov 12007ft . 
allowing us to determine the high-accuracy ICRF positions 
of the target sources through relative VLBI astrometry. 
There is a small overlap (3 common sources) between the 
two fields (Field A around J1616+3621 and Field B around 
J1623+3909). By comparing the images and positions 
derived independently for these common targets, we can 
assess the reliability of our survey method. 

The list of target sources is given in Table [1] for both 
fields. The source names were derived from the VLBI co- 
ordinates wherever available, otherwise the SDSS coordi- 
nates were used. For reference, we list the SDSS spectro- 
scopic redshift (if available) , the SDSS r magnitude and the 
FIRST integrated flux density at 1.4 GHz. We also give the 
5-GHz total flux densities from the GB6 catalogue (Gregory 
et al. I1996|) . Together with the FIRST values, these are use- 
ful to infer the radio spectral slope of the sources. However, 
due to their detection limit of ~ 18 mJy, our weaker tar- 
gets are not found in the GB6 catalogue. The approximate 
positions of the reference and target sources observed are 
shown in Fig. [TJ 
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Table 1. The 13 and 16 sources targeted around the phase-reference calibrators J1616+3621 (Field A, top) and 
J1623+3909 (Field B, bottom) 



Source name 


Equatorial coordinates (J2000) 




z 


SDSS 


FIRST 


GB6 


VLBI 5 GHz 


offset from 






R.A. 


Dec. 








r 


int 


int 


peak 


int 


SDSS position 






h m s 


O 












mjy 


mjy 


mjy/bm 


mjy 


mas 


J161252 


.4+355234 


16 12 52.42647 


35 52 34.3661 






22.0 


21.6 




3.7 


8.0 


23 


J161446 


.9+374607 


16 14 46.95806 


37 46 07.2795 


1 


.53 


16.6 


51.1 


36 


76.9 


85.2 


108 


J161453 


.1+364908 


16 14 53.10423 


36 49 


08.5616 






22.8 


85.8 


35 


6.8 


28.6 


347 


J161550 


.9+364827 


16 15 50.93939 


36 


48 


27.3383 






21.1 


24.3 


22 


4.8 


8.7 


64 


J161606 


.0+363614 


16 16 06.01644 


36 36 14.4327 


1 


.36 


19.9 


24.9 


22 


14.0 


17.9 


33 


J161734 


.3+371154 


16 17 34.34556 


37 11 54.9837 






21.0 


32.6 


- 


4.1 


12.9 


94 


J161748 


.4+380141* 


16 17 4841146 


38 01 41.8062 


1 


.61 


19.2 


21.2 


60 


61.4 


70.4 


53 


J161900 


.7+381139* 


16 19 00.73053 


38 11 


39.3010 






21.4 


44.9 


- 


3.8 


6.5 


51 


J162004 


.7+351554 


16 20 04.73836 


35 15 54.4634 


2 


.94 


18.6 


50.0 


38 


27.5 


38.5 


174 


J162027 


.0+371026 


16 20 27.03237 


37 10 26.8544 






20.4 


77.2 


21 


1.3 


1.4 


156 


J162111 


.2+374604* 


16 21 11.27611 


37 4 6 04.8423 


1 


.27 


18.8 


643.1 


201 


13.5 


46.7 


80 


J162330 


.5+355933 


16 23 30.53332 


35 59 33.1456 





.87 


18.3 


266.8 


144 


83.5 


131.2 


110 


J162419 


.9+353845 


16 24 19.95778 


35 38 45.2701 


1 


.77 


18.2 


28.8 


22 


16.3 


20.0 


47 


J161432 


.1+394445 


16 14 32.19838 


39 


44 


45.0770 






20.6 


316.1 


134 


0.9 


1.6 


675 


J161748 


.4+380141* 


16 17 48.41133 


38 


01 


41.8048 


1 


.61 


19.2 


21.2 


60 


63.2 


70.1 


54 


J161845 


.2+383807 


16 18 45.28153 


38 


38 


07.7917 






19.6 


90.7 


46 


4.2 


8.1 


48 


J161900 


.7+381139* 


16 19 00.73060 


38 


11 


39.3000 






21.4 


44.9 


— 


3.9 


5.6 


50 


J162000 


.4+404319 


16 20 00.47630 


40 


43 


19.1665 






22.4 


103.9 


27 


4.4 


11.9 


107 


J162111 


.2+374604* 


16 21 11.27601 


37 


46 


04.8406 


1 


.27 


18.8 


643.1 


201 


12.9 


20.6 


82 


J162229 


.3+400643 


16 22 29.31193 


40 


06 


43.5915 





.69 


18.6 


32.1 


58 


75.1 


83.5 


88 


J162240 


.7+381637 


16 22 40.72642 


38 


16 


37.3210 






22.5 


73.5 


83 


35.8 


39.8 


48 


J162422 


.0+392440 


16 24 22.00133 


39 


24 


40.9144 


1 


.12 


17.9 


144.5 


128 


87.3 


97.6 


59 


J162453 


.4+375806 


16 24 53.47942 


37 


58 


06.6459 


3 


.38 


18.5 


56.4 


22 


8.6 


23.6 


11 


J162605 


.1+404806 


16 26 05.15018 


40 


48 


06.6508 






22.0 


125.6 




1.7 


5.2 


773 


J162624 


.9+401945 


16 26 24.95879 


40 


19 


45.4248 






19.2 


104.6 


39 


0.7 


0.8 


66 


J162656 


.5+381000 


16 26 56.52929 


38 


10 


00.9853 






20.3 


54.2 


29 


20.1 


24.8 


123 


J162945 


.0+374003 


16 29 45.0826 


37 


40 


03.409 






20.9 


33.0 




< 0.6 






J163020 


.7+375656 


16 30 20.7733 


37 


56 


56.412 





.39 


17.9 


21.7 




< 0.7 






J163302 


.1+392427 


16 33 02.10492 


39 


24 


27.3907 


1 


.02 


16.6 


54.4 


35 


13.1 


19.6 


53 



Notes: Col. 1 - source name derived from the VLBI coordinates (if detected) or the SDSS coordinates; sources marked with * are 



observed in both fields; for these common sources the Field B coordinates are more accurate, but we also list the Field A 
coordinates typed in italics for comparison; Col. 2 - J2000 right ascension (h m s); Col. 3 - J2000 declination (° ' "); Col. 4 - 
spectroscopic redshift; Col. 5 - SDSS r magnitude; Col. 6 - FIRST integrated flux density at 1.4 GHz (mjy); Col. 7 - GB6 total 
flux density at 5 GHz (mjy); Col. 8 - VLBI peak brightness at 5 GHz (mjy/beam), an upper limit is given for non-detections; 
Col. 9 - VLBI integrated flux density at 5 GHz (mjy); Col. 10 - VLBI peak-SDSS position difference (mas) 



3. Observations 

The 5-GHz VLBI observations involving 8 antennas of the 
European VLBI Network (EVN) took place on 2 March 
2007. The total observing time was 6 hours. The partic- 
ipating radio telescopes were Effelsberg (Germany), the 
phased array of the 14-element Westerbork Synthesis Radio 
Telescope (The Netherlands), Jodrell Bank Mk2 (UK), 
Noto (Italy), Toruh (Poland), Onsala (Sweden), Sheshan 
and Nanshan (P.R. China). The recording data rate was 1 
Gbit s -1 . Eight intermediate frequency channels (IFs) were 
used in both left and right circular polarisation. The total 
bandwidth was 128 MHz in both polarisations. 

We observed a calibrator and two programme sources 
in ~ 5-min phase-referencing cycles which included the an- 
tenna slewing and the time spent on the calibrator source. 
Each programme source was observed in 4 scans of ~ 2 min 
duration. These scans were distributed over the whole ob- 
serving period to obtain a better (u, w)-coverage of which 
an example is shown in Fig. [2] The typical restoring beam 
using natural weighting is ~ 4.1 mas x 1.4 mas in a posi- 
tion angle (PA) of —5°. The theoretical image thermal noise 
(la) was estimated as ~ 70 /iJy/beam. The correlation took 
place at the EVN Data Processor at the Joint Institute for 



VLBI in Europe (JIVE) in Dwingeloo (The Netherlands). 
The integration time was 2 s, the undistorted field of view 
was limited to 5" by bandwidth-smearing. 



4. Data processing 

The NRAO Astronomical Image Processing System (AIPS) 
was used for the data calibration (e.g. Diamond 1995). 
Visibility amplitudes were calibrated using system temper- 
atures measured regularly at the antennas during the obser- 
vations. Fringe-fitting was done separately for each phase- 
reference calibrator (J1616+3621, J1623+3909) and other 
secondary calibrator or fringe-finder sources (J1613+3412, 
J1642+3948). The calibrated visibility data were exported 
to the Caltech Difmap package (Shepherd et al. 1994) for 
imaging. The conventional hybrid mapping procedure in- 
volving several iterations of CLEANing and phase (then 
amplitude) self-calibration resulted in the images shown in 
Figs.[3]and[4]for the two phase-reference calibrators. Overall 
antenna gain correction factors were determined and ap- 
plied to the visibility amplitudes of the target sources in 
AIPS. In any case, these corrections were not larger than 
8%. 
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Fig. 1. Sketch of the celestial positions of the two phase- 
reference sources (©) and the 26 programme sources. 
The borders of Fields A and B projected onto the right 
ascension-declination plane arc indicated with dashed lines. 



The fringe-fitting was repeated in AIPS, taking the cal- 
ibrator clean component models into account. The resid- 
ual amplitude and phase corrections resulted from the non- 
pointlike structure of the phase-reference calibrators were 
considered this way. The phase, delay and delay-rate solu- 
tions obtained were interpolated and applied to the target 
source data. The visibility data of the programme sources, 
unaveraged in time and frequency, were also exported to 
Difmap for imaging. Consistently with our analysis in Paper 
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Fig. 2. The (u, w)-coverage of one of the target sources 
(J161252.4+355234) in our 5-GHz EVN observations. 
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Fig. 3. Naturally weighted 5-GHz VLBI image of the ref- 
erence source in Field A (J1616+3621). The lowest con- 
tours are drawn at ±0.5 mJy/beam. The positive contour 
levels increase by a factor of 2. The peak brightness is 
81.4 mJy/beam. The restoring beam is 4.1 mas x 1.5 mas 
at PA=— 7?4. The coordinates are relative to the bright- 
ness peak. The Gaussian restoring beam is indicated with 
an ellipse in the lower-left corner. 



I, for the sources with integrated flux density less than 
15 mjy, the phase- referenced images displayed in Fig. are 
the results of a CLEANing process only. For the brighter 
sources, phase self-calibration with a solution interval of 
30 s was also performed using the CLEAN component 
model. Amplitude self-calibration was only done for the 
sources brighter than ~ 40 mjy. In the images of Fig. [5j the 
coordinates shown are relative to the brightness peak posi- 
tion. However, the technique of phase-referencing preserves 
the sky position of the target sources. The absolute coor- 
dinates of the brightness peaks are given in Table Q] These 
positions were located using the AIPS task MAXFIT. The 
image parameters are listed in Table [2] 

Independently from the phase-referenced solutions, di- 
rect fringe-fitting for 8 strong target sources with peak 
brightness > 20 mJy/beam (Col. 8 in Table [l} was done in 
AIPS. This allowed us to check the images in order to quali- 
tatively verify the source structures. The image parameters 
(peak brightness, noise level) could also be compared with 
those of the phase-referenced images. The dynamic range 
of the latters may be affected by residual errors introduced 
by the switching between the target and reference sources. 
By comparing our phase-referenced images with those ob- 
tained after direct fringe-fitting, we estimated the coheren- 
cence loss. The differences in peak brightness were ranging 
from 4% to 23%, with an average value of 12%. The 5cr up- 
per limits for our two non-detections (J162945. 0+374003 
and J163020. 7+375656) were increased by 25% to account 
for the loss of coherence (Table [TJ . 
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Fig. 4. Naturally weighted 5-GHz VLBI image of the ref- 
erence source in Field B (J1623+3909). The lowest con- 
tours are drawn at ±0.3 mJy/beam. The positive contour 
levels increase by a factor of 2. The peak brightness is 
85.4 mJy/beam. The restoring beam is 4.0 mas x 1.4 mas 
at PA=-3?3. 



5. Results and discussion 

The 5-GHz VLBI images of 24 sources are displayed in 
Fig. [5[ Most of them show nearly unresolved mas-scale ra- 
dio structures, although typical extended one-sided core-jet 
sources are also seen. The sources detected can be grouped 
according to their peak brightness as follows: 

— 5 bright objects (~ 60 — 90 mJy/beam) 

— 7 objects with ~ 10 — 35 mJy/beam 

— 9 objects with ~ 2 — 8 mJy/beam 

— 3 weak objects (~ 1 mJy/beam) 

In comparison with the results obtained in the DEVOS 
pilot experiment (Paper I), not just the detection rate is 
more favourable here (85% vs. 40% of the initial FIRST- 
based sample), but also the ratio of the bright objects 
within the detected sample is considerably higher. 

We note that only about a half of the sources - those 
with r Si 20 m - in our present sample have measured spec- 
troscopic redshifts in SDSS. This is in fact due to the spec- 
troscopic target selection method applied for quasars in 
SDSS (Richards et al. 12002)) . The relation between the opti- 
cal r magnitude (taken from the SDSS DR6) and the 5-GHz 
flux density in the VLBI-detected components is shown 
in Fig. [6] Sources from Paper I that satisfy our present 
selection criteria are also included in the plot. Although 
our sample is rather small for a rigorous statistical eval- 
uation of this relation, it can be seen that the optically 
brighter sources tend to have higher radio flux densities in 
their VLBI components. In particular, the mean 5-GHz flux 
density is 53 mJy for the 12 quasars with known redshift 
(r < 20 m ), while for the remaining 15 quasars (r > 19 m ) 
this value is 23 mJy. Most of our compact sources that show 
nearly unresolved VLBI structures with no extended fea- 
tures (5 out of 6) are drawn from the brighter sub-sample. 




20 19 
SDSS r magnitude 



Fig. 6. The 5-GHz VLBI flux density (mJy) versus the 
SDSS r magnitude of the quasars detected with VLBI (this 
paper and Paper I). Spectroscopic redshifts are available 
for the objects marked with squares. 



The 5-GHz radio luminosity of the 12 VLBI-detected 
DEVOS quasars with known redshift (Paper I and this 
work) is displayed in Fig. [71 The luminosity values arc 
around 10 26 W Hz -1 . Indeed these are one or two orders 
of magnitude smaller than the typical luminosities in VLBI 
surveys of bright AGNs. To illustrate this, we also plotted 
in the same figure the VLBI Space Observatory Programme 
(VSOP) 5-GHz pre-la unch VLB A survey data (VLB Apis; 
Fomalont et al. [2000) for comparison. Here we used the 
sum of the fitted model component flux densities to calcu- 
late luminosities of the sources with redshift values available 
(the overwhelming majority of the VLB Apis sample). Note 
that the VLBApls sample also contains AGNs that are op- 
tically identified with galaxies. The quasar luminosities in 
our small sample of 12 sources are similar within a factor of 
~ 14, despite the broad range of redshifts (0.69 < z < 3.38) 
covered. Qualitatively, there is no obvious difference be- 
tween the compact mas-scale radio structures of our low 
and high redshift quasars. 

5.1. Astrometric accuracy 

The ICRF positional accuracy for the target sources de- 
pends on the phase-reference calibrator source position ac- 
curacy, the target -calibrator angular separation, the angu- 
lar resolution of the interferometer array and the signal- 
to-noise ratio (see also the discussion in Paper I). We con- 
servatively assume a relatively weak source (dirty image 
peak/rms ratio 5:1) at the edge of the field at 2° sepa- 
ration from the calibrator. Since a priori calibrator posi- 
tional accuracies are different for our two fields (1.53 mas 
in right ascension, 1.03 mas in declination for J1616+3621; 
0.37 mas in right ascension, 0.56 mas in declination for 
J1623+3909; Petrov I2"0t77j) . we give the following overall 
values. In Field A, the absolute positonal accuracy of the 
target sources is at least 1.6 mas in right ascension and 
1.1 mas in declination. In Field B, these accuracies are 
0.5 mas and 0.7 mas, respectively. Indeed, the compari- 
son between the coordinates derived independently for the 
3 target sources (J161748.4+380141, J161900.7+381139, 
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Fig. 7. The integrated 5-GHz radio luminosity (W Hz -1 ) 
of the VLBI components as a function of the source redshift 
for the DEVOS sources with known redshift (squares) and 
the generally brighter radio-loud AGNs from Fomalont et 
al. pOOO)) (asterisks). 

J162111. 2+374604; Table Q} located in the overlapping re- 
gion of Fields A and B (Fig. Q]) shows that the diferences 
are within 2 mas in right ascension and within 1.7 mas in 
declination for each source. Due to the above reasons, for 
these common objects we consider the coordinates obtained 
in Field B more accurate. 

The SDSS optical positions for all but 3 of the 24 VLBI- 
detected sources are within 180 mas from the peak positions 
measured in the VLBI images (Col. 10 in Tabic [I]). This is 
in good agreement with what is expected and found for the 
SDSS offsets with respect to the VLBI positions for a large 
sample of common objects (e.g. Frey et al. I2006al [2006b)) . 
The extremely large (~ 0'/7 - 0'.'8) VLBI-SDSS posi- 
tional differences seen in the case of J161432. 1+394445 and 
J 162605. 1+404806 might be due to a significant displace- 
ment between the radio and optical brightness peaks, most 
likely because of an arcsecond-scale extended radio emis- 
sion. Our VLBI positions are in fact closer to (within 0"5 
of) the FIRST radio positions which are accurate to ~ 1" 
only. Another possibility is that we see random associa- 
tions. The chance of finding pairs of unrelated sources in 
our small 26-element sample that are projected so close 
to each other on the sky is not negligible: Ivezic et al. 
(2002) estimate 3% probability for random matches be- 
tween SDSS and FIRST objects within our search radius of 
1'.'5. Therefore we choose a conservative approach and do 
not consider J161432.1+394445 and J162605. 1+404806 suc- 
cessful VLBI detections of the corresponding SDSS quasars. 

5.2. Notes on individual sources 

Our two phase- reference calibrators (Figs. Oil]), as well as 
three of our brightest target sources (J161748. 4+380141, 
J162240. 7+381637 and J162330. 5+355933; Figs. [5k, q and 
1, respectively) are found in the VLBA Imaging and 
Polarimetry Survey (VIPS; Helmboldt et al. [2TT"T7|FI image 
data base. This latter survey concentrates on the celestial 
area covered by SDSS and provides polarization-sensitive 

5 http : //www. phys .unm. edu/~gbtaylor /VIPS 



5-GHz imaging data obtained using the NRAO Very Long 
Baseline Array (VLBA). Their sample of more than 1100 
sources was drawn from the VLA Cosmic Lens All-Sky 
Survey (CLASS; Myers et al. [2003]) . selected on the ba- 
sis of flat radio spectrum and a lower limit of 85 mJy flux 
density at 8.4 GHz. The brightness distributions seen in 
our images are consistent with the VIPS results for all the 
common sources. 

J161446. 9+374607 The compact structure observed 
is consistent with the blazar classification (Marcha et al. 
2001). The 5-GHz flux density we detected is more than 
twice as high as the total flux density found in the GB6 
catalogue. This source (Fig. [5]b) might therefore be highly 
variable. 

J162111. 2+374604 Marecki et al. (f2006|) describe this 
quasar (known also as 4C +37.46) as an asymmetric core- 
jet-lobe source. They imaged it with MERLIN at 5 GHz, 
and with VLBI at 1.7 and 5 GHz as part of a subarcsecond- 
size Compact Steep Spectrum (CSS) radio source sample. 
Our image in Fig.[5js shows the inner jet structure extended 
to the north-east. Guided by the lower resolution images 
of Marecki et al. (2QQ6|), we were also able to marginally 
detect a small fraction (~ 1 mJy/beam) of the diffuse radio 
emission originating from the lobe at around 368 mas east 
and 376 mas north of the core (out of the field displayed in 
Fig. [5k). 

J162229.3+400643 The 5-GHz flux density we mea- 
sured is significantly higher than the GB6 value, possibly 
indicating variability of the source (Fig. [5p) . Neumann et 
al. (|1994j) give S5 ghz = 58 mJy and a flat radio spectrum. 

J162422.0+392440 Two other compact FIRST radio 
sources are found within 13" of this object. These three 
form a nearly linear structure with an angular extent of ~ 
21", the central object being the brightest. This is a doublc- 
lobed FR-II radio quasar (de Vries et al. 2006) of which we 
imaged the core with VLBI (Fig. [SJ') . It appears compact 
with a hint on an extension to the N-NE which might have 
been resolved with longer north-south baselines. 

J162453.4+375806 A high-redshift (z = 3.38) quasar 
with a resolved mas-scale structure in the NE-SW direction 
(Fig. [Sis). Benn et al. (2005) describe this source as the most 
radio-luminous broad absorption line (BAL) quasar known, 
with unusual radio and optical properties. Its radio emission 
is strongly polarised (11% at 10 GHz) and the rotation mea- 
sure (18350 rad m~ 2 ) is among the highest values known for 
any extragalactic source. The GHz-peaked radio spectrum 
of the source suggests a young age. Based on their 22-GHz 
VLA observations, Benn et al. p705)) give 0'.'4±0'.'1 for the 
angular size of J162453. 4+375806. On the other hand, the 
spectral turnover frequency (~ 2 GHz in the rest frame) in- 
dicates ~ 0.1 — 1 kpc linear size. The integrated flux density 
in our VLBI components (23.6 mJy) is practically equal to 
the GB6 value (22 mJy) and the total flux density mea- 
sured by Benn et al. (2005) (23.3 mJy). Thus we can con- 
clude that the entire radio emission comes from a region as 
small as ~ 10 mas, corresponding to a projected linear size 
of less than 0.1 kpc (assuming a flat cosmological model 
with tt m = 0.3, n A = 0.7 and H = 70 km s^ 1 Mpc" 1 ; this 
model is used throughout the paper). 

J163020.7+375656 We did not detect this object 
with VLBI at 5 GHz (peak brightness < 0.7 mJy/beam). 
McLeod & McLeod (|200ip studied the host galaxy of 
this optically luminous low-redshift (z — 0.39) quasar 
which they call radio-quiet with the Near-Infrared Camera 
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and Multi-Object Spectrometer (NICMOS) on the Hubble 
Space Telescope (HST). 

J163302. 1+392427 Another FIRST source is found 
within 9" to the south-east of our object. Its direction is 
close to the position angle of the mas-scale extension seen 
in our VLBI image (Fig. [5^). Intra-night optical variability 
was detected in this quasar at one epoch by Stalin et al. 

pro . 

6. Conclusions 

Phase-referenced 5-GHz VLBI images of 24 radio- 
loud AGNs at mJy flux density level or higher 
are presented here. To our knowledge, 20 of these 
sources have never been imaged with VLBI be- 
fore. Five of our target sources (J161446. 9+374607, 
J161748.4+380141, J162229. 3+400643, J162330. 5+355933, 
J162422. 0+392440) are nearly as bright and compact at 
5 GHz as our phase-reference calibrator sources. It is in- 
teresting to note that two of them (J161446. 9+374607 and 
J162422. 0+392440) have so far escaped from being selected 
for any of the VLBI surveys. This indicates that the ra- 
dio sky is still not fully explored even for relatively strong 
(~ 80 mJy) mas-scale compact extragalactic sources. 

We detected 22 out of 26 target SDSS quasars, 85% of 
our initial sample. (Two other VLBI detections are prob- 
ably chance positional coincidences with SDSS quasars.) 
Combining this with the DEVOS pilot project results 
(Paper I) and counting also the phase-reference calibrator 
sources, the VLBI detection rate is nearly 90% for the ra- 
dio quasars that are known to be compact in both optical 
(SDSS) and in the FIRST survey, and have integrated flux 
density Sxa > 20 mJy in FIRST. Based on currently avail- 
able optical and low-resolution radio sky surveys, we found 
a highly efficient way to pre-select extragalactic sources 
that are detectable with VLBI at 5 GHz. On the other 
hand, for any given object found in the SDSS and FIRST 
lists, the selection method we used is also an excellent tool 
for assessing whether it can be potentially successfully im- 
aged with VLBI. 

The SDSS DR6 contains nearly 9000 quasars that sat- 
isfy our selection criteria and thus can be targeted with 
VLBI. Extrapolating our detection rate for this entire sam- 
ple, the number of potential VLBI detections can be as 
high as 7500. Such a big VLBI-optical sample - partly with 
spectroscopic redshifts available - would eventually be in- 
valuable for a multitude of astrophysical and cosmological 
studies. 

Astrometric projects like global or local densifications of 
the ICRF, and linking the radio reference frame with the 
future optical reference frame would also benefit from the 
DEVOS data. In practice, the data base is to be built up 
field by field gradually, concentrating on celestial areas that 
are otherwise targeted by different projects, e.g. deep spec- 
troscopy, or deep imaging at any waveband. The immediate 
science return of these observations could be enhanced this 
way. 
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Table 2. VLBI image parameters for Field A {top) and 
Field B {bottom) 



Source 


Peak 


Lowest 


Restoring beam 


name 




contour 


size 




PA 




(mjy/beam) 


(mas x mas) 


n 

V / 


J161252. 4+355234 


3.7 


±0.33 


4.2 


X 


1.5 


-8.6 


J161446.9+374607 


76.9 


±0.50 


4.1 


X 


1.4 


-7.4 


J161453.1+364908 


6.8 


±0.39 


4.2 


X 


1.5 


-8.5 


J161550.9+364827 


4.8 


±0.31 


4.2 


X 


1.5 


-7.7 


J161606.0+363614 


14.0 


±0.33 


4.2 


X 


1.5 


-7.6 


J161734.3+371154 


4.1 


±0.42 


4.1 


X 


1.5 


-7.0 


J161748.4+380141 


61.4 


±0.32 


4.0 


X 


1.5 


-9.0 


J161900.7+381139 


3.8 


±0.37 


3.9 


X 


1.5 


-8.1 


J162004.7+351554 


27.5 


±0.30 


4.3 


X 


1.4 


-6.3 


J162027.0+371026 


1.3 


±0.30 


4.2 


X 


1.4 


-5.8 


J162111. 2+374604 


13.5 


±0.30 


3.7 


X 


1.0 


-5.0 


J162330.5+355933 


83.5 


±0.75 


4.1 


X 


1.4 


-5.2 


J162419.9+353845 


16.3 


±0.42 


4.2 


X 


1.4 


-3.9 


J161432. 1+394445 


0.9 


±0.37 


4.2 


X 


2.4 


-25.0 


J161845.2+383807 


4.2 


±0.42 


4.1 


X 


1.4 


-1.2 


J162000.4+404319 


4.4 


±0.39 


4.0 


X 


1.4 


-4.5 


J162229.3+400643 


75.1 


±0.29 


4.1 


X 


1.4 


-2.3 


J162240.7+381637 


35.8 


±0.25 


4.2 


X 


1.4 


-0.2 


J162422. 0+392440 


87.3 


±0.37 


4.1 


X 


1.4 


-1.4 


J162453.4+375806 


8.6 


±0.34 


4.1 


X 


1.3 


-0.8 


J162605. 1+404806 


1.7 


±0.46 


4.5 


X 


2.3 


-18.5 


J162624.9+401945 


0.7 


±0.31 


4.1 


X 


1.3 


-0.1 


J162656.5+381000 


20.1 


±0.38 


4.1 


X 


1.3 


-0.5 


J163302. 1+392427 


13.1 


±0.41 


4.0 


X 


1.3 


3.5 



Notes: Col. 1 - source name; Col. 2 - peak brightness at 5 GHz 
(mjy/beam); Col. 3 - lowest contour level (mjy/beam) corre- 
sponding to 3cr image noise; Col. 4 - Gaussian restoring beam 
size (mas X mas); Col. 5 - restoring beam major axis position 
angle (°) measured from north through east. 
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Fig. 5. Naturally weighted 5-GHz VLBI images of 24 programme sources. The image parameters (peak brightness, lowest 
contour level corresponding to 3er image noise, Gaussian restoring beam size and orientation) are listed and explained in 
Table [2] The positive contour levels increase by a factor of 2. The coordinates are relative to the bightness peak for which 
the absolute coordinates are given in Table [T] The restoring beam is indicated with an ellipse in the lower-left corner. 



S. Frey et al.: Deep Extragalactic VLBI-Optical Survey (DEVOS) II., Online Material p 4 



(f) 



J161734.3+371154 



(g) 



J161748.4+38014J 1 



Right Ascension (mas) 



20 -20 

Right Ascension (mas) 



1 

(h) 


1 1 1 

C 


O 


1 

(i) 


1 1 1 



11 


\ o 













1 ' - 


Declinatioi 













<D 
O 

m 

o 




If 

J 

a 


1 

1 


J161900.7+381139 

I , I 


CM 

1 


1 

1 




J162004.7+351554 





Right Ascension (mas) 



-20 
Right Ascension (mas) 



1 

(j) 


1 1 1 




1 


• 

1 


J162027.0+371026 



(k) 








J162111.2+374604 



Right Ascension (mas) 



-20 
Right Ascension (mas) 



Fig. 5. (continued) 
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Fig. 5. (continued) 
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Fig. 5. (continued) 
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Fig. 5. (continued) 



